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Abstract 

We study the dynamical description of scalar-tensor gravity by performing the best-fit analysis 
for two cases of exponential and power-law form of the potential and scalar field function coupled 
to the curvature. The models are then tested against observational data. The results show that 
in both scenarios the Universe undergoes an acceleration expansion period and the geometrical 
equivalent of dark energy is associated with a time-dependent equation of state. 
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1. INTRODUCTION 



Recent cosmological observations of high redshift type la supernovae, the galaxy clusters 
survey Sloan digital sky survey (SDSS) [5j and Chandra X-ray observatory [6] 

discover the universe accelerating expansion. Also Cosmic Microwave Background (CMB) 
anisotropics observations [7J indicate universe flatness and that the total energy density 
is close to one [8j. The observations determine basic cosmological parameters with high 
precisions and also strongly indicate that the universe presently is dominated by a smoothly 
distributed and slowly varying dark energy (DE) component. A dynamical equation of 
state ( EoS) parameter that is related directly to the evolution of the energy density in 
an expanding universe can be taken as an appropriate parameter to interpret the universe 
acceleration p]-fTO]. 

Suggested by string theories [llj that, the scalar-tensor models provide the simplest 
model-independent description of uniflcation theories which predict couplings between scalar 
flelds and curvature. They have assumed a prominent role in cosmology since any uniflca- 
tion scheme, such as supergravity in the weak energy limit, or inflationary models such as 
chaotic inflation, seem to be supported by them [12]. In addition, they derive the universe 
acceleration and possible phantom crossing p!3]-[T9]. 

In this paper, we study the dynamics of the scalar-tensor theories by flrst best fltting the 
model with the observational data for distance modulus using statistical test. This allows 
us to observationally verify model before checking the results against any other experimental 
data. The layout of the paper is as follows. Section two recalls the basic dynamics of scalar- 
tensor theories. In Section three, in terms of new dynamical variables and for two power 
law and exponential forms of the functions F((/>) and ^(0) in the model, we numerically 
solve the equations. In section four, we best flt the model parameters and initial conditions 
against observational data for distance modulus. In section flve we examine the model with 
two cosmological tests; the Eos parameter and the observational data for Hubble parameter. 
Finally we summarize and conclude in section six. 
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2. THE MODEL 



A general action in four dimensions, where gravity is nonminimally coupled to a scalar 
field 0, is given by 

where R is the scalar curvature of 5^^^^, g its determinant and and ¥{(/)) are two generic 
functions representing the coupling of the scalar field with geometry and its potential energy 
density respectively. We chose WnG = 1. The dynamics of the real scalar field depends 
a priori on the functions F{(j)) and V {(/)). By using the transformation (j) = F{(f) and 
^(^) ^ 2dF((^)/dip action (1), the Brans-Dicke theory can be recovered. Also, the 

standard Newton coupling is recovered in the limit ^ — |. 

The action [T] with negative sign of the kinetic term seems to contain ghost fields for 
certain values of the non-minimal coupling functions F {(/)). The presence of such fields 
would explain how the dominant energy condition is violated in the model, as is required to 
give rise to phantom behavior. Ghost fields in the theory lead to instability of the vacuum 
and the model may suffers from catastrophic theoretical instabilities. However, as follows by 
assuming power law or exponential forms of and best fitting the model, the function 

is always positive as also required for the conformal transformation to be well defined. 
The field equations can be derived by varying the action ([T]) with respect to Qj^jy 

1 

T 

where 



FmR,^--g,^R) = Tl,t\ (2) 



T^t^ = - \9,..M;a - 9^.aF{(l>) + F(0);^, + g^,V{(l>). (3) 

In addition, variation with respect to the scalar field (/) gives the klein-Gordan equation, 

dF 

□0 + + n^) = 0. (4) 

Assuming a spatially homogeneous scalar field (f>{t), and in background of a flat FRW metric, 
the field equations become, 

3H F = -3HF+^ + -, (5) 

jo T/ 

-2HF = 3H^F + 2HF + F+ — + -, (6) 

4 2 

# + 3H4)^ = 6{H + 2H^)F - V, (7) 



where equation ^ is the energy constraint corresponding to the (0,0)-Einstein equation. 
In comparison with the standard cosmological model, one can find an effective EoS for the 
model as p^ff = uJeffPeff^ where Pe// and p^ff are respectively equivalent to the right hand 
side of the equation ^ and ([g]). In general, the system of nonlinear second order differential 
equations ([5])-([7]) has no known analytic solution. However, in the following by introducing 
new dimensionless dynamical variables, we replace the above equations with a system of 
nonlinear first order differential equations which is capable to solve numerically. 



3. DYNAMICAL SYSTEM 



In this section, we study the structure of the dynamical system by introducing the fol- 
lowing dimensionless variables, 

F (j)^ V 

In two different scenarios for the scalar field function F^cp) and potential 1^(0) we investigate 
the dynamics of the model. The cosmological models with such functions have been known 
lead to interesting physics in a variety of context, ranging from existence of accelerated 
expansions [20j to cosmological scaling solutions [2T]-[22]. 
Power law and V{(j)) 

Using equations ([sJ-Q, and power low forms for F{(/)) = (j)^ and V{(/)) = 0^, the dynam- 
ical equations of the new variables become, 

_A = tj^^Q 2VL:—, (10) 

where N = ln{a) and also 

3n.3 + i^4-^. (13) 
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By using the Friedmann constraint equation ^ which now becomes, 

Qv + ^F + ^i = l- (14) 



the equations ([9|)-(11) reduce to two differential equations for and Q^. However, the 
differential equations are long and tedius and so are not presented in here. 
Exponential and 1^(0) 

Similarly, using equations ([sJ-Q, and exponential forms of F{(/)) = e^^ and V^((/>) = e^^ 
the evolution equations become, 

dQi ad 2^1 H 



where 



^ - 3fif + — + - (19) 



Together with the Friedmann constraint equation (14) for this scenario similar to the previ- 



ous case the equations (15)-(17) reduce to two differential equations for and Q^. Again, 
since the differential equations are long and messy, we omit them. 

Still, solving analytically the set of non-linear first-order differential equations in both 
cases of exponential or power-law and ¥{(/)) is difficult. So in the following we present a 
numerical solution for them. In addition, we simultaneously best fit the model parameters a, 
/3 and initial conditions $^^(0), ^^^^(0), H{0) with the observational data using the method 
when solving the equations for two exponential and power-law scenarios. The advantage of 
simultaneously solving the system of equations and best fitting is that the solutions are 
physically meaningful and observationally favored. 

4. BEST-FITTING AND COSMOLOGICAL CONSTRAINTS 

In the following we numerically solve the field equations and also best-fit our model 
parameters with three sets of observational data; i) The Sne la dataset, ii) Cosmic Microwave 
Background (CMB) data, iii) Baryon Acoustic Oscillations (BAO) data. 



The Sne la dataset 

The difference between the absolute and apparent luminosity of a distance object is given 
by, = 25 + 51og;L0 where the Luminosity distance quantity, dL{z) is given by 

By numerical techniques, we solve the system of dynamical equations for and in both 
power law and exponential cases. In addition, to best fit the model parameters and initial 
conditions, two auxiliary equations for the luminosity distance and the hubble parameter 
are employed: 

(d (22) 

To best fit the model for the parameters a, 13 and the initial conditions ^^^(0), ^^^(0), H{0) 
with the most recent observational data, the Type la supernovea (SNe la), we employe the 

method. We constrain the parameters including the initial conditions by minimizing the 

function given as 

1=1 



where the sum is over the SNe la data. In relation (23), jJf^^ and jif^^ are the distance 
modulus parameters obtained from our model and observation, respectively, and a is the 
estimated error of the jif^^ . 
CMB data 

For CMB data, the CMB shift parameter R, given by 



1 /'^^rec 

^^"-i W) '''' 

where E[Z) = is the redshift of recombination z^ec = 1090[25j. The parameter R ties 
up the angular diameter distance to the last scattering surface, the comoving size of the 
sound horizon at z = 1091.3 and the angular scale of the first acoustic peak in CMB power 
spectrum of temperature fiuctuations [23] . [24] . The updated value of R from WMAP5 is 
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Robs = 1.710 ± 0.019 [25]. It should be noted that, apart from the SNIa and BAO data, the 
R parameter can provide information about the universe at very high redshift. The Xcmb 
for the CMB data is 



= (^H^ (25) 

where the corresponding la errors is aR = 0.019. 

BAO data 

For BAO data, from the measurement of the BAO peak in the distribution of SDSS 
luminous red galaxies, we define parameter A as [26j 



where Zb = 0.35. The SDSS BAO measurement [26j gives Aobs = 0.469(n,/0.98). 0.350.017, 
where the scalar spectral index is taken to be Ug = 0.960as measured by WMAP5 [25j. The 
parameter A is nearly model-independent and imposes robust constraint as complement to 
SNIa data. The for the BAO data is 



xIao = ^^^-^ (27) 

where the corresponding la errors is = 0.017. Table I shows the best fitted model 
parameters and initial conditions in both power law and exponential cases. 



TABLE I: Best-fitted model parameters 



Cosmological test 


model 


a P 


npio) H{o) 




SNIea 


power law 


-1.25 4.44 


1.6 


0.6 0.694 


544.6282552 


SNIea + CMB + BAO 


power law 


-1.5 5.95 


1.6 


0.6 0.720 


683.5751181 


SNIea 


exponential 


-1.92 5.88 


1.6 


0.6 0.698 


542.9071824 


SNIea + CMB + BAO 


exponential 


-1.68 6.7 


1.6 


0.6 0.722 


675.4777593 
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FIG. 1: The best-fitted one dimension likelihood and confidence level for a, /3 and Hq for expo- 
nential functions 

Figs. 1-4 shows the constraints on the model parameters a, /3 and the initial conditions 
for J^f(O), Q^{0) and H{0) at the 68.3%, 95.4% and 99.7% confidence levels in both cases 
of power law and exponential functions. 

In Fig. 5, the the numerically calculated distance modulus, fi{z), in the model is best 
fitted with the observational data for the model parameters a, f3 and initial conditions 
for Qf{0)^ ^(/)(0) ^iid H{0) using method in both cases of power law and exponential 
functions. 



Next, we examine the best fitted model against some other observational data. Unfor- 
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SNIea + CMB + BAO 




FIG. 2: The best-fitted two dimension likelihood and confidence level for a, /3 and for expo- 
nential functions 

tunately, the nature of dark matter and dark energy in the universe is such that direct 
observation are not possible. From current deceleration parameter of the universe we are 
able to determine the approximate value of the current EoS parameter. However to quant i- 
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FIG. 3: The best-fitted two dimension likelihood and confidence level for a, ^ and Hq for power 
law functions 

tatively predict the behavior of the EoS parameter of the universe in the past we need more 
observational data. 

A direct verification of our model parameters is a comparison of the best fitted Hubble 
parameter obtained from numerical solution of the field equations with the observational 
data. In the following section, we test our model against the observational data for Hubble 
parameter and also the known dynamical behavior of the EoS parameter of the unvierse. 
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SNIea 





FIG. 4: The best-fitted two dimension likelihood and confidence level for a, /3 and i^o for power 
law functions 
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FIG. 6: The best-fitted effective EoS parameter plotted as function of redshift for F{(j)) and V 
left) exponential, right) power-law 

5. COSMOLOGICAL FITTING RESULTS 

The effective EoS parameter in both cases of power law and exponential functions in 
terms of the new dynamical variables respectively are given by 

With the best-fitted model parameters and initial conditions with the observational data 
in both cases, the effective EoS parameters are shown in Fig. 6. 

As can be seen, a common feature in both cases of power law and exponential is that 
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FIG. 7: The best-fitted H{z) plotted as function of redshift for and V left) exponential, right) 
power-law 




FIG. 8: The best-fitted reconstructed F(0) and plotted as function of redshift for left) 

exponential forms, right) power- law forms 

the best fitted EoS parameters in high redshifts show radiation dominated era. However, in 
power-law case, phantom crossing never occurs in the past or future whereas in exponential 
case, the best fitted parameter with Sne la exhibits phantom crossing in the future and the 
best fitted EoS parameter with Sne Ia+ CMB + Bao shows phantom crossing in the past. 

A second cosmological test with the observational data for Hubble parameter is shown in 
Fig. 7 in both power paw and exponential cases. From the graphs we observe that in both 
cases, our model is in good agreement with the observational data. 

In addition, the behavior of the best fitted reconstructed functions F((/)) and potential 
V{(/)) versus (/) in both exponential and power law cases are also shown in Fig. 8. 
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6. DISCUSSION AND SUMMARY 



This paper is to present the dynamics of the scalar-tensor gravity theories in cosmology. 
We assume two different scenarios for the scalar field function F{(/)) coupled to the geom- 
etry and the potential in the exponential and power law forms. In both cases, by 
representing the model in terms of the new dynamical variables we simultaneously solve the 
field equations in fiat FRW cosmology and best fit the model parameters and also the initial 
conditions with the observational data for distance modulus using the test. Best fitting 
the model parameters provides physically reliable and observationally verified solutions. 
The best fitted effective EoS parameters, ujeff^ in both scenarios exhibit a late accelerating 
universe. 

In power law case, the universe begins from a decelerating state in high redshift (radiation 
dominated era), transits to an accelerating state at about z ^ 1 in the past (quintessence 
phase with cjg// > ~1) again in near future re-enter the deceleration state. In this 
scenario, the universe never experience phantom crossing and phantom phase. Similarly, in 
the case of exponential function parameterizations, the universe begins from deceleration 
state in the high redshift and enter the acceleration state in the near past at about the same 
redshift. However, it will also transit to the phantom phase (cjg// < — 1) in the future for the 
best fitted effective EoS parameter with Sne la data and in the past for the best fitted with 
Sne Ia+ CMB + Bao data. However, in power law case, constrained by only Sne la, the 
EoS parameter in future transits from phantom to quintessence phase and finally to matter 
dominated phase. From the graph we also observe that the current values for effective EoS 
parameter in both exponential and power law cases are between —0.4 and —1.1. 

The phantom crossing feature is impossible to reproduce in the context of single field 
quintessence or even phantom models. It is interesting that in our scalar tensor model with 
phantom field (negative kinetic term), and best fitted power law F(0), the crossing occurs 
while the requirement of F{(/)) > is guaranteed. In this paper, we have also constructed 
experimentally viable function F{(j)) and potential V{(j)) with the best fit parameterizations 
obtained from the recent SNe la, CMB and Bao dataset. . From Fig. 7, the model shows 
a good agreement with the observational data for Hubble parameter. Fig. 8 shows the 
behaviour of the the best fitted F((/)) and V{(j)) for both exponential and power. The best 
fitted reconstructed functions show monotonic increasing/ decreasing behavior in each model 
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with respect to the scalar field 0. The positive F{(f)) in both cases is in agreement with the 
theoretical requirement. 
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